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Abstract 

Our previous study of the near-threshold pp — > ppTr^ reaction based on 
a hybrid nuclear effective field theory is further elaborated by examining the 
momentum dependence of the relevant transition operators. We show that the 
two-pion exchange diagrams give much larger contributions than the one-pion 
exchange diagram, even though the former is of higher order in the Weinberg 
counting scheme. The relation between our results and an alternative counting 
scheme, the momentum counting scheme, is also discussed. 



In the standard nuclear physics approach (SNPA), a nuclear reaction amplitude 
is calculated with the use of the transition operator derived from a phenomenological 
Lagrangian and nuclear wave functions generated by a high-precision phenomenolog- 
ical NN potential. SNPA has been enormously successful in explaining a vast range 
of nuclear phenomena. Meanwhile, a nuclear chiral perturbation approach based on 
heavy-baryon chiral perturbation theory (HBxPT) is gaining ground as a powerful 
tool for addressing issues that cannot be readily settled in SNPA. HBxPT is a low- 
energy effective field theory of QCD, based on a systematic expansion in terms of the 
expansion parameter e = Q/A-^ <C 1, where Q is a typical energy-momentum involved 
in a process under study or the pion mass m.^, and the chiral scale A^ ~ 47r/,r — 1 
GeV. HB^PT has been applied with great success to low-energy processes includ- 
ing e.g., pion-nucleon scattering and electroweak reactions on a nucleon and in few- 
nucleon systems. Our present work is concerned with a HB^PT study of the near- 
threshold pp—yppiT^ reaction. A motivation of this study may be stated in reference 
to the generic NN—>-NN7r processes near threshold. Although HB^PT presupposes 
the small size of its expansion parameter Q/A^, the pion-production reactions involve 
somewhat large energy- and three-momentum transfers even at threshold. Therefore 
the application of HB^PT to the NN^NNir reactions may involve some delicate 
aspects, but this also means that these processes may serve as a good test case for 
probing the limit of applicability of HB^PT. Apart from this general issue to be 
investigated, a specific aspect of the pp—^ppn^ reaction makes its study particularly 
interesting. For most isospin channels, the NN^NNtt amplitude near threshold is 
dominated by the pion rescattering diagram where the vrA^ scattering vertex is given 
by the Weinberg- Tomozawa term, which represents the lowest chiral order contri- 
bution. However, a quantitatively reliable description of the NN^NNtt reactions 
obviously requires detailed examinations of the corrections to this dominant ampli- 
tude. Meanwhile, since the Weinberg- Tomozawa vertex does not contribute to the 
pion-nucleon rescattering diagram for pp^ppn^, this reaction is particularly sensi- 
tive to higher chiral-order contributions and hence its study is expected to provide 
valuable information to guide us in formulating a quantitative description of all the 
NN^NNtt reactions (including the channels that involve a deuteron). 

The first HB^PT-based study of the near-threshold pp^ppir^ reaction was made 
in Refs. [H [2] . In HB^PT one naturally expects a small cross section for this reac- 
tion since, for s-wave pion production, the pion-nucleon vertex in the impulse ap- 
proximation (lA) diagram and the pion-rescattering vertex in the one-pion-exchange 
rescattering (Ivr-Resc) diagram arise from the next-to-leading-order (NLO) chiral la- 
grangian. A remarkable feature found in Refs. [H [2] is that a drastic cancellation 
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between the lA and Ivr-Resc amplitudes leads to the suppression of the pp-^ppn'^ 
amplitude far beyond the above-mentioned naturally expected level. This destruc- 
tive interference is in sharp contrast with the constructive interference reported in 
SNPA-based calculations [31 H]. It is to be recalled that the pp— >pp7r° cross section 
obtained in Refs. [3l H] was significantly smaller (by a factor of ~5) than the experi- 
mental value [3]. The drastic cancellation between the lA and Ivr-Resc terms found 
in the HBxPT calculations [U |2] leads to even more pronounced disagreement be- 
tween theory and experiment. In this connection it is worth noting that, according 
to Lee and Riska [6J, the heavy-meson {a and u) exchanges can strongly enhance 
the pp— >pp7r° amplitude. It is also to be noted that cx-meson-exchange introduced 
in many A^A^ potentials is more properly described by correlated two-pion-exchange 
(see e.g., Refs. [ZIE]), and that there have been substantial developments in deriving 
a two-pion exchange A^A^ potential using HBxPT, see e.g. [9]. These developments 
were conducive to a HBxPT study of two-pion-exchange (TPE) contributions to the 
pp^ppTT^ reaction [TUl [H]. In the plane- wave approximation it was found [10] that 
TPE contributions are indeed very large (as compared to the Ivr-Resc amplitude), a 
result that is in line with the finding in Ref . [H] . A subsequent DWBA calculation [TT] 
indicates that this feature remains essentially unchanged when the initial- and final- 
state interactions are taken into account. More recent investigations p!2l [T3l [T^ [T5]. 
however, have raised a number of important issues that call for further investigations, 
and the purpose of our present note is to address these issues. 

In Ref. [TOj, to be referred to as DKMS, were derived all the transition operators 
for pp— i>pp7r° belonging to next-to-next-to-leading order (NNLO) in the Weinberg 
counting, and these operators were categorized into Types I ~ VII, according to the 
patterns of the corresponding Feynman diagrams; see Figs. 2 - 5 in DKMS. Types I, 
II, III and IV belong to diagrams of the two-pion exchange (TPE) type, while Types 
V, VI and VII arise from diagrams of the vertex correction type. A notable feature 
pointed out in DKMS is that the contributions of Types II ~ IV are by far the largest, 
and that they even exceed those of the Ivr-Resc amplitude, which is formally of lower 
chiral order. On the other hand, the possibility of strong cancellation among the TPE 
diagrams was pointed out in Refs. [121 113] ■ This motivates us to make here a further 
study of the behavior of the TPE diagrams^ 

A remark is in order here on a counting scheme to be used. At the NN^NNtt 
threshold the nucleon three-momentum must change from the initial value p ~ 
^JlrQm^ to zero, entailing a rather large momentum transfer. To take this large 

^For a brief report on this study, see Ref. [I6] . 
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momentum transfer into account, Cohen et al. [2] proposed a new counting scheme, 
to be called the momentum counting scheme (MCS); see Ref. [13] for a detailed 
review. In MCS the expansion parameter is e = p/niN — i'm-,r/^NY^'^, which is 
larger than the usual HBxPT expansion parameter e ~ uiT^/mN- A study based 
on MCS [13] indicates that the Ivr-Resc diagram for pp—>-ppTT^ is higher order in e 
(and hence less important) than a certain class of TPE diagrams, called "leading 
order loop diagrams" , and that MCS is consistent with the estimates of the TPE and 
other diagrams reported in DKMS. Furthermore, according to Hanhart and Kaiser 
(HK) [T2j, the "leading parts" (see below) of these MCS "leading order" diagrams 
exhibit exact cancellation among themselvesH see also Lensky et al. [13]. Although 
these studies are illuminating, we consider it important to examine the behavior of 
the "sub-leading" parts (in MCS counting) of these TPE diagrams in order to see 
whether they can be still as large as indicated by the phenomenological success of the 
Lee-Riska heavy-meson exchange mechanism. In what follows we shall demonstrate 
that this is indeed the case. 

Analytic expressions for the pp—>-pp7r^ transition operators to NNLO in HBxPT 
were given in DKMS. Although these expressions are valid for arbitrary kinematics, 
we find it illuminating to concentrate here on their simplified forms obtained with 
the use of fixed kinematics approximation (FKA), wherein the energies associated 
with particle propagators are "frozen" at their threshold values. In FKA, the TPE 
operator corresponding to each of the above-mentioned Types I ~ IV can be written 
as: 

T = (^^Y^.k)t{p,p',x) (1) 

where p {p') is the relative three-momentum in the initial (final) pp state {pi—p2 = 2p, 
p( — P2 = 2p'), k = p — p', X = p ■ p', and S = |(ai — 0^2). The function t{p,p',x) 
diverges as — 00, and it is useful to decompose t{p,p',x) into terms that have 
definite ^-dependence as — 00. It turns out [18] that t{p,p',x) can be expressed as 

t(p,p',x)'~°° h{gA/i8f^)y\k\+h{\n{\k\VA'}) + ts + 6t{p,p',x), (2) 

where ^3 is asymptotically /c-independent, and 6t{p,p',x) is 0{k~^). For each of 
Types I ~ IV, analytic expressions for tj's (z = 1,2,3) can be extracted [IB] from 
the amplitudes T given in DKMS [10]. The first term with ti in eq.([2]) is the leading 

^HK [12 pointed out that the sign of the contribution of Type II in Ref. [10] should be reversed; 
we have confirmed the necessity of this correction. 
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part in MCS discussed by HK \T2], whereas the remaining terms, which we refer to as 
the "sub-leading" terms, were not considered by HK. The study of these sub-leading 
terms is an important theme in what follows. Table 1 shows the value of ti for Type 
K (K= I ~ IV) extracted from the results given in DKMS. The third row in Table 1 
gives the ratio Rk = Tx/Tnesc, where Tk is the plane- wave matrix element of T in 
eq.([T]) for Type K (K=I ~ IV) normalized by Tfi^sc, the plane- wave matrix element of 
the iTT-Resc diagram. The fourth row in Table 1 gives R^ = T^/TResc, where Tj^ is 
the plane- wave matrix element of T with the ti term in eq.([2]) subtracted. We can see 
from the table that the most divergent ti terms of the TPE diagrams add up to zero, 
confirming the result of Ref . [1^ . However, this does not necessarily mean that the 
TPE diagrams are unimportant, because we still need to examine the contributions 
of the "sub- leading" terms (the t2, H and 5t terms) in eq.([2]). Comparison of Rk 
and R^ indicates that the subtraction of the ti term reduces the magnitude of Tk 
drastically (except for Type I which has no ti term), but the fact that \Rk\ is of the 
order of unity (Types I, II and IV) or larger than 1 (Type HI) suggests that the TPE 
contributions can be quite important. The sum of the contributions of Types I ~ IV 
is 

Y.Rk {=Y.Rk) = -4.65, (3) 

K K 

which indicates that, at least in plane- wave approximation, the TPE contributions 
are more important than the Ivr-Resc contribution. 

Table 1: For the four types of TPE diagrams, K= I, II, HI and IV, the 
second row gives the value of ti defined in eq.(l2]), and the third row gives 
the ratio Rk = TK/Tj^esc, where Tk is the plane- wave matrix element of 
T in eq.([T]) for Type K, and T^esc is the Ivr-Resc amplitude. The last row 
gives Rk = /Tptesc, where is the plane- wave matrix element of T in 
eq.([T]) with the ti term in eq.([2]) subtracted. 



Type of diagrams : K = 


I 


II 


HI 


IV 


{tl)K 





1 


1/2 


-3/2 


Rk 


-.70 


-6.54 


-6.60 


9.19 


Rk 


-.70 


-0.82 


-3.73 


0.61 



Next we investigate the behavior of the TPE diagrams as we go beyond the plane- 
wave approximation by using distorted waves (DW) for the initial- and final-state 
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A^A^ wave functions. For formal consistency we should use the A^A^ potential derived 
from HBxPT, but we adopt here a "hybrid EFT" approach and use phenomenological 
potentials. A conceptual problem in adopting this hybrid approach is that, whereas 
the TPE transition operators derived in HB^PT are valid only for a momentum range 
sufficiently lower than A^~l GeV, a phenomenological A^A^ potential can in principle 
contain any momentum components!! To stay close to the spirit of HBxPT, we 
therefore introduce a Gaussian momentum regulator, exp(— p^/A^), in the initial and 
final distorted wave integrals, suppressing thereby the high momentum components 
of the phenomenological A^A^ potentials; this is similar to the MEEFT method used in 
Ref.[T9]. A(3 should be larger than the characteristic momentum scale of the pp^ppir^ 
reaction, p ~ ^Jmfqm^, ~ 360 MeV/c, but it should not exceed the chiral scale A^; 
in the present study we shall consider the range, 500 MeV< A^ <1 GeV. As high- 
precision phenomenological A^A^ potentials, we consider the Bonn-B potential [20], 
the CD-Bonn potential [21j, and the Nijm93 potential of the Nijmegen group [22]. 

It is worth noting here that several groups [231 Elj have developed a systematic 
approach to construct from a phenomenological potential an effective A^A^ potential, 
called Viow-ki that resides within a model space which only contains momentum 
components below a specified cutoff scale Kiow-k- In this work we will use ^low—k &s 
derived by the Stony Brook group [2l]. It is conceptually natural to use Viow-k in 
conjunction with transition operators derived from HBChPT A problem however 
is yiow-k l24j, primarily meant for describing sub-pion-threshold phenomena, was 
obtained with the use of a rather low cutoff, Kiow-k ~ 2 fm~^. This cutoff is perhaps 
too close to the characteristic momentum scale p ~ 360 MeV/ c for the pion production 
reaction. It therefore seems worthwhile to "rederive" Viow~k employing a momentum 
cut-off higher than 2 fm~^ and use it in the present DWBA calculation. Below we 
will use Viow-k generated from the CD-Bonn potential for Aiow-k= 4 and 5 fm~^. We 
remark that, as is well known, Viow-kS, generated from any realistic phenomenological 
potentials lead to practically equivalent half-off-shell A^A^ K-matrices and hence the 
same NN wave function. 

We evaluate the TPE contributions in DWBA for a typical case of Tiab = 281 
MeV. Since the ti terms in eq.(l2]) add up to zero, we drop the ti terms in our cal- 
culation!! Thus, in eq.(l), we use t*{p,p',x) instead of t{p,p',x), where t*{p,p',x) is 

^ A pragmatic problem associated with this conceptual issue is that, in a momentum-space calcu- 
lation of the matrix elements of the TPE operators sandwiched between distorted pp wave- functions 
generated by a phenomenological NN potential, the convergence of momentum integrations is found 
to be extremely slow [171 HI] • 

^ Removing the ti term lessens the severity of the convergence problem in our momentum inte- 
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obtained from t{p,p', x) by suppressing the ti term. The partial-wave projected form 
of t*{p,p',x) in a DWBA calculation is written as: 



J 



V Stt 



p dpp dp' / dxtpisoip') f{p,p,x) {p-p'x)iJ3p^{p) (4) 



Here ipa{p) is a distorted two-nucleon relative wave function in the a partial- wave 
(^Sq for the initial state and ^Pq for the final state) given by 



^a(p) = COs{Sa) S{p - Pon)/p'^ +'P 



' Ka{p,poj 



(5) 



where 5a is the phase-shift for the a partial wave, and Ka{p,Pon) is the partial- 
wave K-matrix pertaining to the asymptotic on-shell momentum Pon- The plane-wave 
approximation corresponds to the use of the wave functions of the generic form: 



(6) 



We show in Table 2 the values of J, eq.(jl]), for the TPE operators of Types I ~ IV, 
calculated at Tiab = 281 MeV, with the use of the Nijm93 potential of the Nijmegen 
group [22j ^ and Viow-k- For the Nijm93 potential case, we present the results for five 
different values of between 500 and 1000 MeV/c. For the Viow-k case, the results 
for two choices of Kiow-k are shown: Aiow-k = 4 fm~^ and 5 fm~^. For comparison, the 
values of J corresponding to plane- wave approximation are also shown (bottom row). 
From Table 2 we learn the following: (1) The results for the Nijm93 potential with the 
gaussian cutoff are stable against the variation of A^- within a reasonable range 
(500 - 1000 MeV/c); (2) There is semi-quantitative agreement between the results for 
the NijmQS potential and those for Viow-k] (3) A semi-quantitative agreement is also 
seen between the DWBA and PWBA calculations; (4) The feature found in the plane- 
wave approximation that the contributions of the TPE diagrams are more important 
than the Ivr-Resc contribution remains unchanged in the DWBA calculation; the 
summed contribution of the TPE operators is larger (in magnitude) than that of 
Ivr-Resc by a factor of 2~3.5. 



gration mentioned in footnote 3. 

^ We have checked the resuhs obtained using the Bonn-B and CD-Bonn NN potentials are very 
similar to those for the Nijm93 potential case, which we show here as a representative case. 
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Table 2: The values of J, eq.(jl]), corresponding to the TPE diagrams of 
Types I ~ IV, evaluated in a DWBA calculation for Tiab = 281 MeV. The 
column labeled "Sum" gives the combined contributions of Types I ~ IV, 
and the last column gives the value of J for Ivr-Resc. For the Nijm93 
potential case, the results for five different choices of Aq are shown. For the 
case with Viow-k, CD-4 (CD-5) represents Viow-k generated from the CD- 
Bonn potential with a momentum cut-off Aiow~k = 4fm~^ (5fm~^). The 
last row gives the results obtained in plane-wave approximation. 





I 


II 


III 


IV 


Sum 


Ivr-Resc 


VNijm 


Ag = 500 MeV/c 


-0.11 


-0.12 


-0.55 


0.08 


-0.70 


0.18 


H^ijm 


Ac = 600 MeV/c 


-0.12 


-0.12 


-0.57 


0.07 


-0.74 


0.20 


H^ijm 


Ag = 700 MeV/c 


-0.12 


-0.11 


-0.57 


0.06 


-0.74 


0.21 


H^ijm 


Ag = 800 MeV/c 


-0.12 


-0.11 


-0.55 


0.04 


-0.74 


0.22 


^ijm 


Ag = lOOOMeV/c 


-0.12 


-0.10 


-0.52 


0.03 


-0.71 


0.23 


Vion^^k (CD-A) 


-0.12 


-0.09 


-0.46 


0.03 


-0.65 


0.23 


Vion^-k (CD-5) 


-0.09 


-0.06 


-0.30 


-0.01 


-0.46 


0.22 


Plane— wave 


-0.06 


-0.07 


-0.30 


0.05 


-0.37 


0.080 



We now discuss the above results in the context of MCS . A subtlety in MCS is 
that a loop diagram of a given order u in i not only contains a contribution of order u 
("leading part") but, in principle, can also involve contributions of higher orders in e 
( "sub-leading part" ) due to the non-analytic functions generated by the loop integral. 
As mentioned, however, HK [12] considered only the leading part, which correspond 
to the ti term in eq.([2]). According to MCS, for the reaction pp — > ppn^, the loop 
diagrams corresponding to our Type II, III and IV diagrams belong to NLO in the 
e parameter, whereas those corresponding to Type I and the Ivr-Resc tree diagram 
are next order in e (NNLO); see Table 11 in Ref. [13]. Meanwhile, as discussed 
earlier, the sum of the "leading parts" of the NLO diagrams vanishes, and therefore, 
in calculating J's in Table 2, we have dropped the ti term contribution, retaining only 
the "sub-leading" parts of these NLO diagrams. This means that all the entries in 
Table 2 represent "sub-leading contributions" (NNLO) in MCS. If we look at Table 2 
from this perspective, we note that the order-of-magnitude behavior of our numerical 
results is in rough agreement with MCS, although Type IV tends to be rather visibly 
smaller (in magnitude) than the others. However, it is striking that J for Type III 
is significantly (if not by an order of magnitude) larger than the other sub-leading 
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contributions. (A similar feature was also seen in R* in Table 1.) In view of the fact 
that Type III arises from crossed-box TPE diagrams [TU], there is a possibility that 
the enhancement of the Type III diagrams may be related to the strong attractive 
scalar A^A^ potential that is known to arise from TPE crossed-box-diagrams [71 18] . 

We have studied the "sub-leading" parts, which are of NNLO in the momentum 
counting scheme (MCS) [T3], of the TPE amplitudes for the pp^ppir^ reaction in both 
PWBA and DWBA calculations. We have shown in fixed kinematics approximation 
(FKA) that, even though the leading parts of the TPE amplitudes cancel among 
themselves fT^i [H], the contributions of the sub-leading parts are quite significant. 
They are in general comparable to the Itt-Rcsc amplitude, and the sub-leading part 
of the Type III diagrams is even significantly larger than the Ivr-Resc diagram. The 
total contribution of the TPE diagrams is larger (in magnitude) than that of the 
Ivr-Resc diagram by a factor of ~5 (PWBA) or 2~3 (DWBA). We have focused here 
on the TPE loop diagrams but, to obtain theoretical cross section for pp^ppn^ that 
can be directly compared with the experimental value, we must consider the other 
diagrams discussed in DKMS as well as the relevant counter terms. These will be 
discussed in a forthcoming article [2^ . 

The authors are indebted to Christoph Hanhart and Anders Gardestig for useful 
discussions. A helpful communication from Ulf Meissner is also grateful acknowl- 
edged. This work is supported in part by the US National Science Foundation, Grant 
No. PHY-0457014, and by the Japan Society for the Promotion of Science, Grant-in- 
Aid for Scientific Research (C) No. 15540275 and Grant-in- Aid for Scientific Research 
on Priority Areas (MEXT), No. 18042003. 
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